Defects, order, and hysteresis in driven charge-density waves 



0^ 



o 

O 

o 



Mikko Karttunen^, Mikko Haataja^, K. R. Elder^, and Martin Grant-^ 
^Department of Physics and Centre for the Physics of Materials, McGill University, 3600 rue University, Montreal (Quebec), 

Canada H3A 2T8 
^Department of Physics, Oakland University, Rochester, MI, 48309-4487 

(February 1, 2008) 

We model driven two-dimensional charge-density waves in random media via a modified Swift- 
Hohenberg equation, which includes both amplitude and phase fluctuations of the condensate. As 
the driving force is increased, we find that the defect density first increases and then decreases. 
Furthermore, we find switching phenomena, due to the formation of channels of dislocations. These 
results are in qualitative accord with recent dynamical x-ray scattering experiments by Ringland et 
al. and transport experiments by Lemay et al. 
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The effect of quenched disorder on a periodic medium 
has been the subject of intense study during the past 
decade. The rich phenomena displayed by these systems 
has raised intriguing questions [Q-||. These include the 
nature of the depinning transition, the effect of disorder 
on structural properties, and the possibility of the ex- 
istence of nonequilibrium analogs of solids and liquids. 
Examples of such systems include charge-density waves 
(CDW) Ql, flux line lattices in type-II superconductors 
lH,|j, and magnetic bubble arrays Q. A common effect 
of quenched disorder in equilibrium systems displaying 
periodicity is to reduce or destroy that periodicity, by 
pinning the system in a disordered state. However, a 
sufficiently large driving force in a nonequilibrium sys- 
tem can act to depin the system, thereby reducing the 
effect of quenched disorder to that of annealed disorder, 
analogous to thermal noise. 

A CDW is a periodic modulation of the electron den- 
sity, of wavenumber Q, which results from the electron- 
phonon interaction S. The local excess density of the 
conduction electrons is given by p{x, t) = Pc{x, t) cos{Q ■ 
X + 4>{x,t)), where Q = 2fcF, kp is the Fermi wavenum- 
ber, pc is the CDW amplitude, and is the phase. 
It is convenient to introduce a scalar order parameter 
i){x,t) ^ Re[e^Q-^pc{x,t)e"f'^^'*^]. 

An important question concerns the role of ampli- 
tude fluctuations of the condensate close to the depin- 
ning transition H-O]. A common approximation is to 
neglect amplitude fluctuations, leading to a phase-only 
description p4[-p7[|. However, this approximation is not 
always justified, as demonstrated by Coppersmith and 
Millis p|-pd| for a model of a driven CDW, and as ar- 
gued by Balents and Fisher fl^. Moreover, recent theo- 
retical studies in 2D indicate that dislocations proliferate 
in the presence of quenched disorder even in the absence 
of an external driving force |1^. Therefore, the role of 
amplitude fluctuations close to the depinning transition 
warrants a closer examination. 

To this end, we propose a model for driven CDWs, 



and show that dislocations indeed proliferate close to the 
depinning transition. This behavior has profound impli- 
cations on the structure of the CDW and its transport 
properties. In particular, we show how the dynamical 
generation of dislocations can make the system more dis- 
ordered above the depinning threshold. This result is in 
agreement with recent x-ray scattering experiments by 
Ringland et al. ||l[ who observed an increasing correla- 
tion length upon letting the CDW relax from the sliding 
state to the pinned state. This behavior, generic in our 
model, has not been seen in any of the previously studied 
theoretical models pg|-|21[|. Furthermore, our results sug- 
gest that upon further increasing the drive, the system 
becomes more ordered. This is in accord with very recent 
experimental results by Ringland et al. B. In particular, 
we predict that the correlation length should increase ex- 
ponentially. These results are consistent with the picture 
of Balents and Fisher ||l^ of a dynamical phase transition 
into a temporally periodic "moving solid" phase. Finally, 
we show that the behavior of the dislocation density can 
be hysteretic upon increasing and decreasing the drive. 
Indeed, we find that the hysteresis is associated with the 
activated formation of channels of dislocations. Based 
on these results we argue that the current in the model 
should display "switching" (i.e., hysteresis). This is in 
qualitative agreement with recent transport experiments 
by Lemay et al. yj. 

We propose a phenomenological equation of motion, 
related to the Swift- Hohenberg equation [g2[, which per- 
mits periodic solutions to occur only in the x direction, 
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where ip = Tp{x,t) denotes the local condensate density, 
and e = |T — Tpj/Tp, with Tp denoting the critical tem- 
perature of the Peierls transition. E denotes the exter- 
nal drive, and V^ is the localized impurity pinning, qc 
stands for the wavenumber of the periodic modulation 



which is set to unity hereafter. Impurities are character- 
ized by a concentration c and quenched Gaussian corre- 



lations {v^^^ 



and {V^^pVi^p 



) = 2Vb5ij. The use of 



'imp I ^ ^^'.^ \' imp* imp I ^^O'^ij- 

both impurity strength Vq and impurity concentration c 
to describe impurities provides a convenient way to tune 
between weak and strong pinning regimes, as described 
below. We note that the dynamics of ip, in the absence 
of external drive, follow from the Lyapunov functional: 
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That is, dip/dt = —dF/6ip. It is straightforward to in- 
clude the effects of thermal noise [|3| . 

This equation naturally produces periodic structures, 
permitting the study of local fluctuations in the wave- 
length, as well as order and defects, as described by 
the phase and amplitude variations of the charge-density 
wave. In particular, a second-order derivative is suffi- 
cient in the y-direction since no periodic ordering occurs 
in the that direction, while both second and fourth-order 
derivatives are necessary in the x-direction to produce 
a periodic pattern, and to suppress spatial gradients on 
very small length scales. Higher-order derivatives in the 
X 01 y directions may alter the local structure of the 
defects, but will not qualitatively change the long-range 
interaction between defects. In the limit of small e, it 
is equivalent to the standard Ginzburg-Landau approach 
P,0-[l6|. The equivalence can be shown via a multiple- 



scales expansion |23 using the slow variables ei, e 



1/2 



X, 



and e^/^y, and the ansatz, V' = e"'^/^V'o + e'0i + £'^^^V'2-t-- • ■■ 
To leading nontrivial order, we obtain the equation of mo- 
tion for the complex amplitude A for the most unstable 
mode dA/dt = A^id"^ /dx^ + d"^ /dv^)A + A-'i\A\^A, which 
is the time-dependent Ginzburg-Landau equation of mo- 
tion for a complex order parameter, which includes both 
phase and amplitude fluctuations. The effects of driving 
force and pinning potential can be straightforwardly in- 
cluded. This, together with an explicit mapping of our 
approach to that of others, will be given in a future pa- 
per P3| . The advantage of our approach is the natural 
interplay between local variations in wavelength during 
the dynamical process, which is difflcult to isolate in the 
Ginzburg-Landau equation. Away from threshold, e ~ 1, 
other unstable modes can become important. We expect 
that, since our equation captures the correct symmetries 
and essential physics of the driven charge density waves, 
that it will remain a valid description for large e. Indeed, 
Qc is only weakly dependent on e [M. Also, our numeri- 
cal results and tests for varying driving forces and e show 
no qualitative dependence on the size of e. 

The nature of CDW pinning is determined by compar- 
ing the energy gain the wave experiences by adjusting 
to the local impurities, to the increased elastic energy 
due that adjustment [p5|. The former energy is Vbce-^/^, 



where e^/'^ gives the amplitude of the CDW. The latter 
is ec£~^, where e gives the elastic force constant, and 



-l/d 



is the distance between impurities in d di- 



mensions. In two dimensions, the ratio of these gives a 
dimensionless parameter k = Vo/(ce^'^). Strong pinning 
corresponds to large k, weak pinning to small n. Herein 
we conveniently tune between regimes by varying c, so 
that strong pinning corresponds to small c. 
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FIG. 1. (a)-(d): Typical configurations. In (a), the exter- 
nal drive E ^ Q, Vo = 0.4 and c = 0.8. In (b), S = 0.05. 
In (c), E = 0.1. In (d), E = 0.2. Only one quarter of the 
original system size is shown for clarity. 

We studied the discretized version of Eq. (|l|) on a 
d — 2 square lattice of size L = 256, where the dis- 
cretized elements of time and space were dt = 0.05 and 
dx — 0.78, respectively. Gradients were evaluated using 
standard finite differencing. We verified that making dt 
and dx smaller did not change our results. The control 
parameter, e, was varied between 0.1 and 1.2, the pinning 
strength Vq — 0.1 — 1.5, the concentration of pinning sites 
c = — 1, and the driving force E — Q — 0.2. Results 
involve averages over at least 10 independent representa- 
tions of the disorder. Typical configurations of the sys- 
tem are shown in Fig. 1 with c — 0.8 and Vq = 0.4. In the 
relaxed state {E = 0), Fig. la, some defects are visible. 
As the driving force is increased past the depinning tran- 
sition (Fig. lb with E — 0.05) the number of defects nr, 
increases |2^. However, for a still larger driving force, 
defects begin to disappear (Fig. Ic with E = 0.1) and, 
eventually, the system becomes more ordered (Fig. Id 
where E = 0.2). 

We quantify this observation by examining the spher- 
ically averaged structure factor S{q), defined by S{q) = 
(|V'(g')P)|g-/|=g, where -0 denotes the Fourier transform 
of tp and (. . .)\qi\=q denotes an average over all configu- 
rations and Fourier modes with \q'\ = q. A typical struc- 
ture factor is shown in Fig. 2. for a relaxed state i? = 0, 
and for increasing driving force in the sliding state where 
E > Ec {E ^ 0.025 and 0.05). As the drive is turned 
off, the peak of the structure factor increases and, cor- 



respondingly, the half width at half maximum decreases. 
This implies that the correlation length increases upon 
switching off the drive, in qualitative agreement with re- 
cent dynamical x-ray scattering experiments by Ringland 
et al. lit]. Furthermore, upon monotonically increasing 
the drive we observe a sharpening of the peak, due to 
decreasing no- This surprising nonmonotonic behavior 
is in good qualitative accord with recent x-ray scattering 
experiments by Ringland et al. 
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FIG. 2. Structure factor S{q) for, from left to right, E ^ 0, 
E = 0.075 and E = 0.09. Vo = 0.5 and c = 0.8. 



relation length exponents, respectively. Indeed, our data 
is in good qualitative agreement with these arguments, 
although our data do not permit a detailed test. 
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FIG. 3. Number of dislocations no as a function of the 
driving force E for strong pinning. Pinning strength Vo ~ 0.5 
and impurity concentration c — 0.4. Circles: increasing 
drive. Squares: decreasing drive. Inset: nn for weak pin- 
ning (c = 0.8). 



To further quantify this nonmonotonic dependence, we 
estimated numerically the number of dislocations njj in 
the system in the following manner. First we determined 
the number of areas ua in the system where the local 
transverse gradient |VyVI > 0.35. We have checked that 
riA = A ■ no + B, where ^ « 10 is the contribution from 
an isolated defect, and i3 is a measure of the roughness of 
the CDW in the absence of defects, which is dependent 
on c. In Fig. 3 we show the number of dislocations njj for 
c = 0.4, corresponding to the strong pinning regime; the 
weak pinning result (c = 0.8) is shown in the inset. As 
Fig. 1 indicated, there is an initial increase in nu^ up to 
a maximum appearing close to the depinning transition, 
followed by a subsequent decay for E > Ec- 

The behavior of n^ and the structure factor can be 
understood using a simple argument. In the absence of 
thermal fluctuations, the generation of dislocations must 
depend on the number density of pinned sites Up and on 
the local strain rate. In addition, annihilation of disloca- 
tions depends on the mobility D and no- In mean- field 
theory, the rates for generation and annihilation of de- 
fects are given by Rq ^rip-v and Ra ~ D ■ nn . The ve- 
locity V of the CDW enters into Rg since the local strain 
rate is proportional to the local velocity. Using the above 
expressions, the mean-field steady-state number of dislo- 
cations is simply no '-^ rip ■ v. Scaling arguments imply 
^ no ^ [E - EcY-" as E ^ Ec and no '- exp-i^^ 
a,s E ^ Ec, where 9 and i' denote the velocity and cor- 
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FIG. 4. The behavior of the system near the depinning 
threshold in the strong pinning regime. Each of the snapshots 
is a difference between two configurations 500 time steps dt 
apart (25 time units since dt — 0.05). The black areas indicate 
where the changes occur, i.e., which areas are moving The 
grey scale is the same in all the pictures. Vo = 0.5, c — 0.4. 
(a) a.t E = 0.020, (b) E = 0.025, and (c) E = 0.030. 

In the strong pinning regime, additional behavior is 
evident: njj shows significant hysteresis as the drive is 
first increased and then decreased. In the weak pinning 
limit, no significant hysteresis is seen. The origin of this 
is evident in Fig. 4. We compare two configurations of the 
same system at different time steps, around the depinning 
threshold. The configurations were taken 500 time steps 
dt apart (25 time units since dt = 0.05), and then they 
were deducted from each other. The darker areas indicate 
where the changes are the largest. Below the threshold, 
in the strong pinning regime, there is little movement 
in the system. As the driving force is increased to the 
transition, a conducting channel appears. As the driving 



force is increased further, the system depins coherently. 
In other words, depinning in the strong pinning regime 
is mediated by activated channels p7[ |. 

In contrast, in the weak pinning regime (not shown) 
domains form cooperatively, even before depinning takes 
place, and the transition takes place collectively, with 
insignificant apparent activation. Unfortunately, in our 
two-dimensional simulations, we cannot determine if the 
qualitative transition in behavior between weak and 
strong pinning is sharp, or governed by crossover. The 
qualitative difference is easy to understand. In the 
strong-pinning regime, impurities are far apart, so chan- 
nels, where regions of current flow, "nucleate" between 
impurities. In the weak pinning regime, where there are 
more impurities, there is insufficient room for channels to 
form, and, instead, a cooperative and collective effort of 
the entire system is required to depin the system. This 
bistable behavior in the strong pinning regime will man- 
ifest itself in the current-voltage characteristics of the 
sample. This is known as "switching", and is a generic 
feature of low-temperature CDWs [U||-|3^. This behav- 
ior, which follows naturally from our approach, is in qual- 
itative agreement with recent transport measurements of 
Lemay et al. |2|. 

In conclusion, we have introduced a novel approach 
to study the of driven charge-density waves. We have 
demonstrated that, above the global depinning thresh- 
old, phase slips appear and lead to an even more disor- 
dered structure for the CDW, as seen experimentally [Q . 
However, far away from threshold the system becomes 
more ordered, as observed very recently 0. These latter 
results are in agreement with the picture of Balents and 
Fisher p^. Finally our model gives rise to a simple and 
generic mechanism for switching in the strong pinning 
regime, through the activated appearance of channels, 
which qualitatively explains recent experimental results 
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